Defect-Suppressed Atomic Crystals in an Optical Lattice 
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We present a coherent filtering scheme which dramatically reduces the site occupation number 
defects for atoms in an optical lattice, by transferring a chosen number of atoms to a different internal 
state via adiabatic passage. With the addition of superlattices it is possible to engineer states with 
a specific number of atoms per site (atomic crystals), which are required for quantum computation 
and the realisation of models from condensed matter physics, including doping and spatial patterns. 
The same techniques can be used to measure two-body spatial correlation functions. 
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There has been a vast amount of recent interest in the 
study of Bose-Einstein Condensates (BECs) Jlj and de- 
generate Fermi gases in optical lattices 0, LalU ■ Such 
systems have many potential applications in quantum 
computing, and also make possible the study of strongly 
correlated systems from condensed matter theory with 
unprecedented control over system parameters. 

However, the study of strongly correlated systems or 
entangled atoms in optical lattices requires not only that 
the corresponding (Hubbard) Hamiltonian 0,0 be prop- 
erly engineered, but also that the system is placed in 
a well defined (pure) initial state with high precision. 
Whilst is possible to load a BEC into an optical lattice 
in the Mott Insulator (MI) phase, which corresponds to 
having a definite number of atoms at each lattice site 
y, l| , non- ideal conditions will always result in defects 
in that phase (i.e., missing atoms and overloaded sites). 
The removal (or precise control) of such defects is then a 
necessary condition for the application of atoms in optical 
lattices to quantum computing and the study of strongly 
correlated systems. In this letter, we propose a coherent 
filtering scheme which, beginning with an uncertain num- 
ber of atoms in each site, provides a method to transfer 
a definite number of atoms at each site into a different 
internal state, and so load a new lattice of atoms with an 
exact number of particles per site. This process dramat- 
ically reduces the site occupation number defects, and 
can be extended, under experimentally reasonable condi- 
tions, to allow the production of doped or pattern loaded 
initial states with almost unit fidelity (i.e., atomic crys- 
tals) . The scheme is also applicable to the production of 
high fidelity initial states in fermion systems (including 
doped states and the loading of composite objects), and 
when combined with moving optical lattices, can be used 
to measure two-body correlation functions. We illustrate 
this by presenting an example of how this scheme could 
be used to both produce and characterise superconduct- 
ing states in an optical lattice with a chosen filling factor. 

We begin by considering a system of bosons loaded into 
an optical lattice such that they do not tunnel between 
neighbouring sites. The atoms are in a particular inter- 



nal state, \a), and have onsite interaction strength (with 
h = k B = 1) U a = ATia s J d 3 x|w(x)| 4 /™, (valid for U a 
less than the separation of the two lowest motional states 
at each site) where a= is the scattering length and w(x) is 
a Wannier function [3| . We then couple the atoms into a 
second internal state, \b), (which is trapped by a second 
lattice potential) via an off-resonant Raman transition 
with Rabi frequency f2(f), which is detuned from state 
| b) by S(t). If we denote the onsite interaction between 
particles in state |6) as J7& and the onsite interaction be- 
tween particles in different internal states as U a b, and 
follow the derivation in we obtain the Hamiltonian 

H = -y-n (n - 1) + -^n&(nt, — 1) + U ab n a n b 

-^-(atb+tfa)-S(t)n b , (1) 

where a and b are the annihilation operators for particles 
in states \a) and \b) respectively, and n a = a^a, fib = wb. 
We can then write the state of the site as \n a , rib), where 
n a and rib are the number of particles in the internal 
states | a) and \b) respectively. We know that the initial 
state at a particular lattice site is \N, 0), where N S 
1,2,..., TVmax. Our goal is to transfer exactly one particle 
from state \a) to state |6), so that for any N, the final 
state is \N — 1,1). In order to do this we must choose 
the initial and final values of the detunings, Si and Sf 
so that the system evolves along the avoided crossing in 
the energy eigenvalues, undergoing an adiabatic passage 
from | JV, 0) — ► \N — 1, 1), and does not evolve along any 
other such avoided crossings. These values have to be 
simultaneously chosen for all values of N, as shown in 
Fig. Once the system is in state \N — 1,1), we can 
turn off the lattice trapping state \a) , leaving a pure state 
with exactly one atom per lattice site. 

The relative locations of the avoided crossings in the 
energy eigenvalues are determined by the values of U a /Ub 
and Uab/Ub- In the adiabatic limit (6 — > 0), we can 
set 51 = 0, and for a given N max it is then straight for- 
ward to find the parameter range for which an appro- 
priate choice of Si and Sf is possible. We consider the 
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FIG. 1: Plots showing avoided crossings in the energy eigen- 
values for N — 1, 2, 3 (Note the vertical scale variation). We 
must choose Si and <5 / (dashed vertical lines) so that we cross 
only from \N, 0) — »• |JV — f , I), as illustrated in the inset. 



locations at which the energy eigenvalues of H cross in 
this limit, and then attempt to find a parameter range 
for which the crossings between \N, 0) and |JV — 1,1) 
for all N fall within a range of detunings [Si , S /] which 
isolates them from all other crossings. The resulting 
boundaries for the allowed parameter space are given by 
Uab = (1 - h)U a + hUb, for i = 1, 2, where for U ab > U a , 
h = 1/AUax, k 2 = l/(2iV max - 4), and for U ab < U a , 
h = l/(2iV max - 2), k 2 = l/(iV max - 3). These bound- 
aries are shown in Fig. [2^ for N max = 4. Larger values 
of iV max result in a more restrictive allowed parameter 
range. Substantially different values of U a and Ub are re- 
quired, which is possible using either spin-dependent lat- 
tices [5j, where U a and Ub are independently controlled 
via the different lattice shape for atoms in each internal 
state, or Feshbach Resonances @, near which the scatter- 
ing length a s is different for atoms in different internal 
states as we tune an external magnetic field. U a b can 
also be independently controlled by slightly displacing 
the lattices trapping internal states \a) and \b). 

In the case of finite-time laser pulses, the adiabatic 
crossings have a finite size, and must be crossed on a 
timescale dictated by the separation of the energy eigen- 
values at the crossings. From the Landau-Zener formula, 
we can estimate the transfer error, e, for \ZN max Q(t) < 
U a ,U b as e w exp[-7riVr2 2 /( 2 cj)]; where 8 ~ \$ f - 8 1 \/t 
and r is the total transfer time. For a large range of pa- 
rameter values, we find that e is very small for r values 
well within the decoherence time (including spontaneous 
emissions, collisional losses, and external parameter fluc- 
tuations, 0) of the optical lattice. This has been con- 
firmed by extensive numerical calculations of the exact 




FIG. 2: a) The parameter range for which appropriate values 
of Si and 5/ can be found when jV max = 4, in the adiabatic 
limit (light shading) and from numerical simulations with a 
smoothed rectangular pulse which rises and falls with a sin 2 (t) 
shape, t — lOOf/fj" 1 and max fi = O.SUb, giving a transfer er- 
ror e < 1% (dark shading), b) Initial and final state fidelities, 
T a and Tb, as a function of temperature, T, for an initial MI 
state described at each site using the GCE. The inset shows 
values of 1 — T a (solid line) and 1 — JF 6 for N a — 2 (dashed) 
and N a = 3 (dotted) on a logarithmic scale. 

time evolution of (Fig. I2K) - 

This filtering scheme naturally cannot correct occu- 
pation number defects in which no atoms are present 
at particular lattice site in state \a). Because in prin- 
ciple extremely low transfer errors, e, are obtainable, 
this is the fundamental limit for the filtering scheme. 
The initial single-site fidelity of a state where we choose 
to have N a particles in each lattice site is defined as 
T a = (N a , 0\w o \N a , 0), where w is the initial density 
operator for the site. Similarly, the final state single- 
site fidelity (for 1 transferred atom) is Tb = (1|u>&|1), 
where Wf, — Tr a w, and w is the final density operator for 
the site. In each case the overall fidelity for M sites is 
given by Tm = T M . If the probability of zero occupa- 
tion in a particular site is po, then Tb = (1 — e)(l — Po)- 
Thus to obtain the highest possible fidelities we should 
start in a MI phase with an average filling factor of 2 
or more particles per site, where po is small. For such 
a MI phase, with an average filling factor of N a at tem- 
perature T, we model the particle number fluctuations 
using a Grand Canonical Ensemble (GCE). For integer 
N a and low initial temperature, exp[— U a /{2T)\ -C 1, we 
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obtain p n — exp[— U a N a /(2T)], so that for e — * 0, Tb ~ 
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(1 — exp[— U a N a / (2T)\). Under the same conditions, we 
can write T a = p(N = N a ) » 1 - 2exp[-Z7 /(2T)]. 
We see (Fig. |2t) that coherent filtering increases the 
fidelity by several orders of magnitude, even when the 
fidelity T a is already reasonably high. Moreover, as 
(1 — Fb) = [(1 — F a )/2] Na , we see that the error in the ini- 
tial state is exponentially suppressed with exponent N a . 
For example, if we assume e — > 0, an initial state with 
N a = 2 and a defect at every tenth lattice site, T a = 0.9, 
T/U a = 0.17, results in a final state with T b » 3 X 10~ 6 , 
i.e., less than one defect every three hundred thousand 
sites. In an attempt to reduce number fluctuations even 
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further we also considered a regime in which the atoms 
were allowed to move between neighbouring sites during 
the filtering process. However, for a fixed mean filling 
factor and transfer time this produces a lower final state 
fidelity than performing the filtering in the MI regime. 

Taking Poisson statistics instead of the GCE, which 
corresponds to the BEC being loaded by suddenly switch- 
ing on the lattice (this can be seen by computing the re- 
duced density operator of the original superfiuid state), 
and assuming a perfect transfer (e = 0), we obtain 
Th = 1 — exp(—N a ). Thus, to achieve high fidelities 
T > 0.99 in the final state requires N a ~ 5. Such N a 
values might ultimately be limited by the increase in 3- 
body collisional losses with many atoms on each site. 

This filtering can clearly be extended to transfer multi- 
ple atoms into state |6) by choosing 8j so that the system 
evolves along more than one avoided crossing (but the 
same number for all initial N) . Similarly, it is possible to 
choose Si and Sf so that the transfer \N, 0) — ► \N — 1, 1) 
occurs only for initial states with a particular value of 
N. Thus, we can measure the defects in the initial state, 
for example by mapping sites with N a = N a + 1 or 
N = N a — 1 atoms onto occupied sites in state \b). 

Spatial patterns can be created by adding a super- 
lattice to the lattice trapping atoms in state \b), site- 
dependently shifting the energy level of \b) and prevent- 
ing atoms in particular sites from coupling to the Raman 
transition. This allows the engineering of many high- 
fidelity spatial states because of the wide range of possi- 
ble superlattice configurations which can be formed using 
overlapping laser beams. 

Coherent filtering can also be applied to a degenerate 
Fermi gas 0,0 in an optical lattice (Fig. J^,). We con- 
sider spin up and spin down fermions (in states \a T, J.)) 
loaded into the motional levels of each lattice site ac- 
cording to the Pauli principle. Atoms in any chosen mo- 
tional level (normally the lowest, as for low temperature 
and high filling factors it will almost always be occu- 
pied) may then be transferred by adiabatic passage to 
a different internal state, \b T)J-)- An example of ap- 
propriate states can be taken from the 2S 1( / 2 level of 
6 Li in a strong magnetic field (i.e., in the Paschen-Back 
regime). We choose \a T, 4) = |m s = 5; mi = 0,-1) and 
\b t,|) = |m s = — i;mi = 1,0), where m s and mi denote 
the magnetic spin quantum number for the electronic 
and nuclear spins respectively. (Where required, spin- 
dependent lattices for these states could be implemented 
for the short times required by coupling the ground states 
between the 2Pi/ 2 and 2P 3 / 2 levels.) As was the case for 
Bosons, we can apply a superlattice to the lattice which 
traps state \b), and hence engineer very high fidelity spa- 
tial patterns (Fig. Ob). By applying different superlat- 
tices during the transfer of different spin states, the ratio 
of spin up to spin down particles in the final state can 
also be very precisely controlled. 

These techniques, when combined with the ability to 




FIG. 3: The production of (c) a BCS state at exactly half 
filling in an optical lattice from (a) a degenerate Fermi gas 
using pattern loading techniques described in the text (a— >b). 



move spin-dependent optical lattices [5j, can be used 
to measure two body correlation functions. Currently, 
sinlge-body correlation functions can be studied using 
the interference patterns produced when atoms are re- 
leased from the lattice, and many-body correlation func- 
tions for a single site can be estimated from collisional 
loss rates. However, for fermion systems in a lattice, 
fourth order correlation functions of the form piji — 
(c| + /+c|_j, Cj_ijCj+;f) (where c%& is the annhilation oper- 
ator for an atom of spin a in site i) are of particular inter- 
est because they provide information about both about 
the long range order and the pair correlation length in 
the system. Specifically, m = Y,i(4+n4~ii^-n^+il) 
characterises the correlation as a function of the sep- 
aration 21 between spin up and spin down particles, 
and 7 fc = Y,t{^\-k^\-ki^+ki^i+k\) gives us information 
about the long range order as a function of the sepa- 
ration 2k between the pairs. These two functions can 
be measured by firstly making the lattice potential very 
deep, and performing the measurement on a timescale 
shorter than that on which atoms can move to neigh- 
bouring sites. Then, to measure Ki, we first transfer all 
of the spin down atoms from their current state, \b I) 
via an adiabatic passage to the internal state \a [)■ Wc 
then shift the lattice trapping the state \a J.) 21 sites to 
the left. If we now transfer the particles in state \a J.) 
to state | b J.) selectively, based on the condition that no 
particle exists (at a particular site) in the state \b t), 
then the number of atoms left in state \a I) is propor- 
tional to the value of m. This is made possible by the 
onsite interaction between particles with different spins 
shifting the detunings at which avoided crossings occur 
in the energy eigenvalues. A similar process can be used 
to measure 7;, but requires an additional step in which 
we flip the spin of particles in state \b [) — > \b f). Whilst 
these processes would be difficult to implement experi- 
mentally, they can be seen as providing an eventual way 
to perform a complete tomography of pair correlations in 
fermion systems. 
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FIG. 4: The pair correlation function pi,j,o, illustrating the 
crossover from MI^BCS states. In (a) (V = 10 J, U = -10 J, 
MI regime) and (b) (V = 0, U = -0.1J, BCS regime) we 
see localized pairs and delocalised pairs respectively. In- 
formation about these correlations can be obtained from 
7i = '^2 i pi_i >i+ i. These results are from numerical diagonali- 
sation of (Sj for 8 sites and half filling with periodic boundary 
conditions. 
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FIG. 5: Numerical values of Ki, which describes the pair cor- 
relation length. These results are examples from numerical 
diagonalisation of (3J for 40 sites and 1 particle of each spin 
type and 16 sites and 2 particles of each spin type (inset) with 
a)V = 0,U= -10J and b) V = 0, U = -0.1J. 



lattices) , consisting of pre- formed pairs in the odd num- 
bered sites (Fig. We slowly decrease V to a small 
value, V < J <C U, allowing the pairs to delocalise. For 
small V but large \U\, the pairs behave as a gas of hard- 
core Bosons. Following the discussion in 0, we find that 
the system is protected throughout this process by an en- 
ergy gap of V. This resulting state is characterised by a 
very short inter-particle correlation length, measureable 
through the correlation function ki (Fig. Et). 

We then lower the onsite interaction until \U\ < J, 
so that the spin up and spin down particles can move 
independently. A meanfield treatment ^3 shows that 
even in the weak coupling limit the state is protected by 
an energy gap which is no less than V. Once the BCS 
state of Cooper pairs sets in Q, we can set V — > by 
completely removing the superlattice. The system thus 
undergoes a transition into a BCS state (Fig. 0t, Et), 
characterised by its long pair correlation length, again 
measureable using m (Fig. |SJd). 

In summary, the coherent filtering scheme presented 
here allows the production of high-fidelity atomic crys- 
tals, and the measurement of important correlation func- 
tions. These features will greatly enhance the application 
of atoms in optical lattices to quantum computation and 
to the precise modeling of condensed matter systems. 
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Such a characterisation would be particularly useful 
in the study of the BEC - BCS crossover in an optical 
lattice ||. Along with initial state preparation, this is 
well illustrated by the method to study the MI to BCS 
state transition in fermions at exactly half-filling shown 
in Fig. |2| We begin by transfering particles of both spin 
types from a degenerate Fermi gas loaded into an opti- 
cal lattice (Fig. via adiabatic passage to a different 
internal state, whilst a superlattice is used to offset even- 
numbered sites in the final state by a potential V (Fig. 
Ob). This system is then described by the Hamiltonian 

N V N 

H = -J C s j> + U^n^fiu + — ^2{-iyn i<a , 

(2) 

where denotes all nearest neighbour combinations 

and hi ta — c icr Ci ;Cr . U denotes the onsite interaction be- 
tween particles of different spin, which we choose to make 
initially large and negative, using a combination of lat- 
tice parameters, and, if necessary, a Feshbach resonance 
. Thus we prepare a state with exactly half-filling (note 
that other filling factors are possible with different super- 
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